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Radiations produced by neutrino-antineutrino annihilation at the pole can be used to heat up 
the primary stage of a thermonuclear warhead and can in principle detonate the device remotely. 
Neutrino-antineutrino annihilation can also be used as a tactical assault weapon to target hideouts 
that are unreachable by conventional means. 
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Nuclear weapon is the most destructive kind among weapons of mass destruction. Hiroshima and Nagasaki are 
lessons in history that shall never be repeated. Since the end of World War II, world leaders had tried to control 
the proliferation of nuclear weapons by political means such as the Nuclear Non-proliferation Treaty in 1968. Many 
04 ' countries did not sign the treaty. In fact it seems that more and more countries are pursuing nuclear weapon programs 
nowadays. After September 11, the concern is that nuclear weapons will fall into the hands of terrorists. Strategically 
speaking the importance of a counter nuclear weapon may soon rival that of the nuclear weapon itself. The purpose of 
f~| , this paper is to explore the possibility of a neutrino counter nuclear weapon technology. The idea of using neutrinos to 
Qh' detonate or melt a nuclear weapon was first proposed by H. Sugawara, H. Hagura and T. Sanami [l|. Their futuristic 
design is based on a 1 PeV neutrino beam operating at 50 GW. It is unlikely that such an intense ultra high energy 
neutrino beam can be realized in the near future. Even if such a neutrino beam is made available, its radiation 
bJO hazard will render it politically nonviable. This paper proposes an alternative idea for a neutrino counter nuclear 
^ weapon that shares some similarities with the previous idea but is technologically feasible, relatively cheap and safe. 
O The present idea is to focus a neutrino beam and an antineutrino beam together in a small region to allow them to 
' ^ , annihilate so that high energy radiations are released as reaction products. The radiations cause neutron spallation in 
the sub-critical nuclear material and initiate fission reactions. The plutonium heats up, ignites the chemical explosive 
i-S^ , around the fissile (fissionable material) in the primary stage of a thermonuclear warhead and subsequently detonates 
the nuclear weapon. The reason of thinking about neutrino for this application is that neutrino cannot be shielded. It 
[ can hit a target such as a nuclear submarine from the other side of the globe and can penetrate a deep underground 
ff^ ■ concrete bunker and missile silo. The trade-off of developing a counter weapon is the introduction of a new weapon. 
^ ' If the new weapon is less destructive than the original weapon, an ethical argument can be made in support of its 
, development. If remote detonation of a nuclear weapon is made possible by a neutrino counter weapon, a nuclear 
■ weapon in the homeland becomes a liability so that there is a real strategic incentive to reduce the stockpile. In that 
: case, there will be a much more convincing political reason to promote non-proliferation. This work aims to study 
• ' the theoretical feasibility of the neutrino counter nuclear weapon as a first step in this direction. The use of neutrino 
as a tactical assault weapon will also be discussed. 

00 : 

\ II. THEORETICAL PREPARATION 

. One of the properties of neutrino is that it has a vanishingly small interaction cross section so that it is the most 
$H ' penetrative radiation known to be in existence today. The invariant amplitude M. is proportional to the weak boson 
. 5t 1 (massive spin- 1) propagator [Hll 



q2 _ + It ' 

where q is the transfer momentum, M is the weak boson mass and the infinitesimal e moves the pole away from the 
branch so that Eq. ([1]) is mathematically well defined. The weak bosons are and . The corresponding masses 
are Mw = 80.403 ± 0.029 GeV and Mz = 91.1876 ± 0.0021 GeV [1]. In the case of an unstable particle with a total 
decay width F, it is customary to make the replacement AP M"^ + iMT according to the Breit-Wigner resonance 
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formula so that Eq. ([T]) becomes @ 

q2 _ 2 _ iMT 
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It is said that the propagator in Eq. ^ violates gauge invariance and that slightly more elaborate modifications are 
needed 0, [!]■ It is also suggested that Eq. ^ is not theoretically justified in quantum mechanics and quantum field 
theory and is motivated mostly by phenomenology Q. Eq. ^ also violates Lorentz invariance in a subtle way. The 
decay width T is typically measured in the rest frame of the particle and is related to the lifetime of the particle tq 
in its rest frame as tq = 1/r. Lifetime is not Lorentz invariant and can be dilated (r = jtq). Traditionally F in 
Eq. ^ is taken to be a constant just like M. As it will be shown later, the width term can be eliminated from the 
final result by special relativistic considerations. For the sake of simplicity, only Eq. ([1]) with e = will be used in 
the calculations of the invariant amplitude Ai in this section. When ^ Af^, the propagator in Eq. ([T]) is reduced 
to ig^v/M'^ and is weighted down by a very heavy weak boson mass such that the invariant amplitude becomes very 
small. It is the reason why neutrino is so non-interactive even though the weak coupling constants gw and gz are 
both larger than the electromagnetic coupling constant ge 01 j 

= iift' 

9z = , (4) 



9e 


sin 6 






9e 


sin 6 


w COS 9 




f 9w 


8 


\Mw 



(5) 

where ge = 7.297352568 x 10"^, 9^, is the Weinberg angle (sin^ 9^^ = 0.23122) and Gp is the Fermi constant {Gp 
1.16637 X 10~^ GeV~^) Q. As a corollary, large Mw is the reason why the inverse beta decay cross section is 
vanishingly small. The present work focuses only on the tree-level diagrams such as vv — > II and vv — > // (Z 
for "lepton", / for "fermion") as shown in Fig. [1] Reactions such as vv — > 77 is possible but is a higher-order 
diagram [l3| . It turns out that a pole is impossible for a charged current in the Feynman diagram of Figure [1^ 
because the transfer momentum is always null {q^ = 0) or spacelike [q^ < 0) as it will be shown later. In the case 
of neutral current in Figure [T|3, the transfer momentum can be timelike (q^ > 0) so that a Z° pole is possible at 
q^ = M|. The Z° pole has been experimentally observed in the spectrum of e+e" scattering at LEP and SLAG since 
1989 Q and in pp scattering experiments at Fermilab [ll|. This work takes advantage of the pole to maximize 
the neutrino-antineutrino annihilation cross section to produce high energy radiation in a maximally efhcient way. 

For the charged current in Figure [1^, The flavor of the outgoing lepton I (antilepton I) must match that of the 
incoming u (i/) at each vertex. Mixed flavors are possible in a reaction as a whole. For example, i^ei^e e^e+, 
i^ei^fi e~ 11'^ and all other combinations arc allowed. The invariant amplitude of II in Fig.[T^ is 

9w \-, „,un „,5\„, 1 I „ '^^9'^ \ \;r,-^M(^ „,5\ 



The transfer momentum is q — pi — p^ = Pp ~ Pi- The approximations made in the second step of Eq. ^ are the 
small masses of neutrino (m^ ^ 0) and lepton (to; <^ Mw)- The invariant amplitude is squared, summed over final 
spins and averaged over initial spins. This calculation is simplified by applying the usual Casimir's trick and trace 
theorems. Neutrino has hclicity h = —1 but antincutrino has h — +1 so that there is only 1 spin state for each initial 
neutrino or antincutrino. The Casimir trick sums over all initial and final spin states by default. By "averaging over 
initial spins" , it simply means that the sum must be divided by the proper factor to avoid over-counting. This way 
the sum over initial spins must be divided by 2 for each initial neutrino or antincutrino. For a reaction involving 
2 incoming neutrinos, the proper invariant amplitude square {Ad"^) needs to be divided by 4. With the standard 
procedure outlined above, Eq. ([5]) can be squared, summed and averaged to give 
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{pi-Pp)iPp-PT)- (7) 



Since the present calculation is focused on high energy (g^ M"^), the rest mass of neutrino and lepton can be 
neglected (m^ — > and mi —f 0). In the lab frame, the 4-momenta can be parametrized as p = {E, p^, Py, Pz) in such 
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a way that 4-momentum is conserved by definition: 

pi ~ £^1(1, sin 6*1, 0, cos 6*1), (8) 

P2 ~ £2(1, -sme2, 0, COS6I2), (9) 

P3 ~ -£3(1, sin0cos0, sinflsini/i, cosfl), (10) 

P4 ^ P1+P2 -~P3- (11) 

The 2;-axis in Eqs. ((8)l-(|lip is chosen to point in the direction of the sum of the initial 3-momenta pi + P2- The 
outgoing 3-momenta P3 and P4 can rigidly rotate around the z-axis through the azimuthal angle 4> and still conserve 3- 
momentum. For the sake of illustration, Fig.[2]shows a special case in that all of the incoming and outgoing 3-momenta 
lie on the same plain. The distance between the foci and JF' corresponds to the vectorial sum Pi + P2 = Pa + P4 
so that 3-momentum is conserved automatically. Assuming the high energy condition Ei ^ rrii, Ei ~ |pi| for 
i € {1, 2, 3, 4}. The ellipsoid £ constrains the point O in such a way that |pi| -I- |p2| — jps] + |p4| so that energy is 
conserved by Ei + E2 — E3 + E4. In essence, the construction in Fig. [2] guarantees the conservation of 4-momentum 
as long as the high energy condition is met. In the special case of 6'i ^ and 62 0, the approximations Ei ~ |pi| 
made in Eqs. ([8l)- (|lip may become invalid and should be replaced by the exact parametrizations 

Px = {Ei, IpI^ sing., cos |p| - sin^; sin;/).,, |p|^cos6'i). (12) 

However small values of 9i and 62 do not have any practical advantage and will not be used in this work. For the 

charged current diagram of Fig. [T^, pi — Pv, P2 — Pv, Pz = Pi and pa — pj- 

The azimuthal angle (j) is not constrained by kinematics and is completely random. The transfer momentum square 
in the charged current diagram of Fig. [T^ can be expressed in terms of the parametrizations in Eqs. ([5|)~ (|lip as 

= {pi - p^Y — -'2E^Ei (1 - sin 9i sin 61 cos - cos 9i cos 6) . (13) 

Regardless of the value of </>, the right-hand side of Eq. (fT3|) is always non-positive. Therefore in (Al^) of Eq. ([7]) 
is either null (g^ — 0) or spacelike {q^ < 0). A pole is impossible in principle in the charged current case because 
it is always true that ^ My^. 

The invariant amplitude for the neutral current diagram in Fig. [TJd is 
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The flavor of uv cannot mix. The type of // can vary but the flavor cannot mix. The values of the neutral vector 
coupling Cy and axial- vector coupling Ca are given by the GWS (Glashow-Weinberg-Salam) Model and are tabulated 
in Table [D Applying the usual rules, the invariant amplitude square can be computed as 



(-Ml) = \ 
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{ {Cl + C\) {{p, ■ psf + (p. • psf\ + ICaCv {{P, ■ Pff - {p. ■ Pff] } . (15) 



The parametrizations in Eqs. ((5|)- (llip can be re-applied with the substitutions pi = Pu, p2 = pp, Ps = Pf and p^ — pj 
The transfer momentum square in the neutral current case is 

g2 = {p,+pi,f ~ 2{p, -p^) = 2E,E^[l ~ cos(0i + 62)]. (16) 

According to Fig.[T)3, < 6*^ < it/2 for i e {1, 2} so that > in Eq. Therefore a pole is possible for (A^|) 

in Eq. (|15p at q^ = . On the average, (coscf)) = and (cos^ 0) = 1/2. The remaining kinematic factors in 
of Eq. (fT^ can be calculated as 



(p^ -Pff = El E'j [1-2 cos 6*1 cos -f - sin^ Oi sin^ 9 -f cos^ Oi cos^ 6 ] , (17) 



{Pi,-Pff = ElE)[l-2cose2Cose+'^siii^92sm^e + cos^92Cos^9y (18) 

Letting Ey = Ej^ leads to 9i = 02 which in turn gives {p„ ■ p/)^ = (pp • pfY so that the form of Eq. (fT5|) is greatly 
simplified. Letting Ei, = Ep has an additional advantage that the intensities of i/g, and Vr in the neutrino beam will 
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match those of Vf. , and Vj- in the antineutrino beam despite of neutrino oscillation as long as the original neutrino- 
antineutrino beams have the same intensities and the matter effect of the earth is cither negligible or canceled by 
symmetry. 

The cross section formula of the reaction l + 2^3 + 4 + -- - + n can be calculated by the Fermi's Golden Rule 0], 



da = {M'' ^ 



{2ny2E3j \{2'Kf2E^J \{2ttY2E, 



x(27r)V(pi+p2-P3-P4 Pn), (19) 

where 5 is a statistical factor that includes a factor of 1/ j for each type of j identical outgoing particles. For the 
Feynman diagrams in Fig. [1] 5* = 1 and n = 4. The differential cross section for the neutral current diagram can be 
obtained from Eq. (fTO|) after performing integrations over pa and P4, 



da _ 1 (M 



dVl (87r)2 E^E„ [1 - cos(6'i + 62)] 

X . (20) 

E^ + Ec,- cos6l^(£;^sine'i - sin 6*2)2 + [E^cosOi cos 6*2)2 

Eq. ([20|) is useful for checking pathology in the theory. In this case, there is no pathology related to the cross section 
formula. It is checked that da/dO, cx Although the exact form of the cross section formula is unimportant 

for the the purpose of this work, Eq. (|^D|) is useful for calculating the relative distributions of the outgoing particles 
once the infinity of l/[q^ — M^] at the pole is canceled out by division. At the pole, the cross section in 
Eq. (j20p is infinite and the interaction probability is 1. This way, the flux of annihilation- induced radiations is 
directly proportional to the intensities of the vv beams. For the purpose of designing an efRcient neutrino weapon, 
Eq. used in = is the key result of this section. The branching ratios of various annihilation-induced 
products // can be calculated by comparing (Al^) using Eq. (|15p . In the case i?^ — Ejy, the branching ratios are 
particularly simple and can be obtained from ratios of Cy + C\ alone. 

Finally the earlier claim that the width term iMT can be eliminated from the final result by special relativity needs 
to be explained. Historically the width is incorporated in the scattering amplitude f{E) in non-relativistic quantum 
mechanics as in [l^] 

It is observed in low energy inelastic scattering experiments that an incoming particle form a compound nucleus with 
the nucleons of the target nucleus as an intermediate state. The discrete energy levels Eq of this compound nucleus 
give rise to resonances. If the compound nucleus is unstable, the resonance is modified as Eo ^ Eo — iT /2. It is 
the so-called Breit-Wigner resonance. In experiments, Eq and F are extracted from partial wave analysis. In hadron 
physics, r is modeled as a function of energy F — > F(W^), where W is the energy in the center-of-momentum (CM) 
frame, along with other model dependent parameters [13|. With minor adjustments, the Breit-Wigner formula is 
phenomenologically robust. Sometimes there are disagreements among different analyses on some of the parameters. 
An example is the i?o+ parameter in the cross section of the 77 photoproduction of 5ii(1535) |3, [isl. fl6l. \vf\ . In 
this case, the discrepancy is thought to have come from different experimental biases and the difficulties in resolving 
nearby resonances. So far the validity of Eq. (PT|) is not questioned. Although the form of Eq. (PT|) has been discussed 
in many quantum mechanics textbooks [T^ , the iMT term in the propagator of Eq. ([5]) in quantum field theory is not 
derived from first principle but is borrowed from quantum mechanics. Reference [6| contains a heuristic explanation 
of Eq. (m as 

(22) 



p2 _ ^4-2 _^ j^/r 2E^{p^ -E^ + i{M/Ep)T/2) ' 

In Eq. ([22|) . is the pole and Ep is the energy of the particle. The factor M/Ep is the reciprocal of the Lorentz 
factor 1/7 that accounts for the relativistic correction of the lifetime of the particle r = jtq. The right-hand side of 
Eq. ([221) is a combination of Eq. (PT|) with the replacement F ^ r/7 and a factor of l/2ii^p that commonly appears 
in relativistic quantum mechanics and quantum field theory. Since the relativistic effect is already included by virtue 
of incorporating I/7 on the right-hand side of Eq. ([22]). F on the left-hand side can be taken as a constant much like 
the rest mass M is a constant. The implicit Lorentz factor 7 on the right hand side of Eq. is frame dependent. 
The left hand side is manifestly frame independent. This ambiguity is a proof that Eq. (I22p is not a theory but a 
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phenomenological patch. Since Eq. (|22p is borrowed from Eq. ([2T|) . it is reasonable to assume that the validity of the 
former is based on the satisfaction of the same conditions of the latter. The basic assumptions of Eq. (PT|) are the 
target nucleus and compound nucleus at rest in the lab frame and a total width F measured in the same rest frame. 
In high energy scattering experiments such as e~^e~ — *■ Z'^ and pp — > Z'^, the equivalent of a target nucleus can be 
taken to be the center-of-mass of the in-state. So far e+e^ and pp collider experiments are always arranged in such a 
way that the beams are colliding head-on with equal and opposite momenta. The CM frame is the natural reference 
frame for the data analysis of such collider experiments. The 3-momentum of Z^ in the CM frame is zero (q = 0) so 
that Z'^ is at rest in this frame. Additionally the center-of-mass frames of the in and out states are also at rest because 
q = 0. Since all the relevant systems are at rest in the CM and lab frames, it is reasonable to expect that the widths 
at the Z^ pole measured by collider experiments agree with the non-relativistic Breit-Wigner formula. In neutrino 
weapon applications, it is seldom practical to shoot the vD beams head-on. In the lab frame where the lyp beams are 
colliding at oblique angles, the momenta of center-of-mass of lyD, Z^ and // are non-zero, q = Pi + P2 = Ps + P4 7^ 0. 
In order to understand the phenomenology of the pole in the q ^ case, the focus is directed to the type of 
collisions in that a resonance is moving relativistically with respect to the target. According to Eq. (fT6|) . the Z^ pole 
is possible in a stationary target when the beam energy reaches 4 x 10^ GeV for a positron beam on an electron target 
or 2 X 10"^ GeV for a anti-proton beam on a proton target. Such experiments do not exist. In the lower energy regime, 
nuclear physics experiments often shoot electron or photon beams at stationary targets. The resonances are usually 
excitations of proton so that the target and the resonance are at rest with respect to each other. In rare high energy 
reactions e+e" Z^Z^ [H, ilSj, ^21,] and pp Z^Z^ [l^l, the condition q ^ is obtained automatically by virtue 
of integrating over p3 and P4 in the calculation of the cross section formula. However Z^ is the out-state in this case 
and not the propagator so that the Breit-Wigner formula does not apply. To the best of the author's knowledge, 
there is no scattering experiment that satisfies the condition of q 7^ at the Z^ pole at the time of the writing of this 
paper. 

At the Z^ pole, most of the energy of the in-state is converted to the Z° mass so that the residual kinetic energy is 
relatively small. The Z^ propagator is approximately at rest in the lab frame even though its momentum is non-zero 
q 7^ 0. Electron and neutrino, on the other hand, have very small mass so that a relatively small amount of transfer 
momentum q will cause them to move relativistically. Z° is a virtual particle. Its influence is limited to a short time 
between the in and out states. As far as the in and out states in their center-of-mass rest frames are concerned, a 
moving Z^ of total width F is effectively the same as a hypothetical stationary of width F/7. The advantage 
of modeling the propagator as a stationary Z° in the center-of-mass rest frame of the in-state is that this condition 
satisfies the assumption of the non-relativistic Breit-Wigner formula in that both the target and resonance are at rest 
with respect to each other. The similar argument can be repeated for the out-state. Following the same reasoning of 
Eq. (|^ , the propagator in the case of q 7^ is taken to be 

(23) 



((72 - Af| + 2imF) ' 

where m e {vriy, my}. The factor of 2 in the width term imF of Eq. (|23p comes from the fact that there are 2 
particles of mass m in the in and out state. The subtlety of Eq. ([22|) is that the argument is crafted first in a preferred 
reference frame (such as the lab frame) and then generalized to be frame independent. Eq ([23|) exploits the same 
strategy by incorporating a Lorentz factor 7 ~ q^/ {2rn) in the center-of-mass rest frame of the in-state (out-state), go 
as in g = (go, q) is canceled from Eq. ([23]) by the factor l/{2qo). With Eq. ([23]) . {M.\) for q 7^ can be re-evaluated 
by averaging the invariant amplitude squares computed with the new total widths corresponding to the in and out 
states, 

{Ml) (X 




MlY + 4m2F2 (g2 _ MlY + 4m2F2 



(24) 

The last step of Eq. (|24|) is motivated by 1/to2 ^ 1/to2. A third term constructed from the rest frame of Z° is not 
needed because it is redundant. In actuality, mi, is small but non-zero. Assuming that - 0.1 eV, M^/mlV'^ - 102^. 
Instead of an infinite cross section as claimed earlier, the reaction rate at the Z^ pole is only 27 orders of magnitude 
larger than that of the inverse beta decay. The improvement is adequate to make the idea of a neutrino counter 
nuclear weapon practical with reasonably intense vi> beams. Therefore the invariant amplitude square in Eq. (|15[) is 
a good approximation up to a factor of 2. 
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III. SIMULATION RESULTS 

The exact details of the designs of thermonuclear weapons are classified. However the basic principles of the 
operation of various types of nuclear weapon are unclassified and freely available through public information [2^ . 
Generally speaking, a modern nuclear weapon is made up of a fission bomb (the primary stage) and a fusion bomb 
(the secondary stage). The primary stage typically has a plutonium core surrounded by a tamper, a neutron reflector 
and a set of chemical explosive lenses on the outermost layer. The tamper slows down the explosion just enough to 
increase the efficiency of the chain reactions. The neutron reflector increases the criticality of the fissile so that the 
efficiency of the explosion is further increased. In some designs (such as the W87 warhead), the tamper and neutron 
refiector seem to be combined in one layer made of beryllium. The secondary stage is a spherical or cylindrical 
complex of solid lithium deuteride (fusion fuel) interlaced with layers of uranium walls and the so-called "spark plug" 
in the center. A spark plug is a fissile lining that enhances compression. The sequence of explosions begins with the 
simultaneous ignition of the chemical explosive lenses in the primary stage to compress the sub-critical plutonium 
core to reach super-criticality by implosion. The heat and pressure generated by the fission bomb in the primary stage 
compress the fusion fuel in the secondary stage. Neutrons from the fission bomb convert the lithium-6 isotopes in 
solid lithium deuteride into tritium (3Li-|-n He-|-f H-l-4.8 MeV). Tritium and deuterium fuse to create a hydrogen 
bomb explosion. In a thermonuclear weapon, most of the energy is generated by fusion in the secondary stage. A 
sketch of the W87 warhead found online [24!| is the source of the geometry used in the simulations of this work. More 
detailed discussions of thermonuclear weapon designs can be found on the same archive f25l|. 

Neutrinos cannot be shielded and can reach a nuclear weapon anywhere on earth or in space as long as the location 
is specified and if the neutrino weapon can focus the high energy vv beams on the target accurately. High energy 
radiations produced by vv annihilation cause neutron spallation inside the fissile material that in turn initiate fission 
reactions. Fissions create heat. When heat is generated in sufficient quantity, it ignites the chemical explosive lenses 
surrounding the plutonium core in the primary stage of the thermonuclear warhead. This way vv annihilation provides 
a means to detonate a nuclear weapon remotely. The critical question is whether these annihilation-induced radiations 
can be generated efficiently and cheaply to make the idea practical. To this end, the simulation tool MCNPX2 .5.0 is 
used to calculate the heat deposition by fissions from high energy external sources at = 45 GeV, approximately 
half the value of Mz- MCNPX2.5.0 is an integration between MCNP4B and LAHET2.8. It uses a very old version of 
FLUKA (FLUKA87) as a hadron generator above the INC (intranuclear cascade) region (E > 10 GeV). The manual of 
MCNPX2 .5.0 explicitly states that it does not include any nuclear reactions for muon. Muon-induced neutron spallation 
has been observed experimentally and is already included in the more recent versions of FLUKA. A newer version of 
FLUKA is planned to be integrated in the next release of MCNPX so that this problem will hopefully be fixed in due time. 
Another problem of MCNPX2.5.0 inherited from LAHET is that not all particle types can be used as primaries. So far 
it is found that MCNPX will crash when t is used as a primary. During the course of this work, it is also discovered 
that MCNPX has a strange bug when 7r° is used as a primary. The bug is already reported to the MCNPX team. In 
principle, spallation neutrons induced by fj,^, r and 7r° can be simulated by the latest version of FLUKA and then 
input into MCNPX to calculate the fission deposition energy. Unfortunately the license agreement of FLUKA restricts 
the software to non-weapon-related use. Since the purpose of the present work is merely a feasibility study and not a 
detailed weapon design, the absence of , t and tt" as primaries in the simulations will not affect the conclusion. The 
remaining particle types used as primaries are e^, tt^, K^, Kl and Ks- Of course neutrinos can also be annihilation 
products. But they will not be a useful source of weapon related radiations and will not be simulated. Despite of 
the shortcomings aforementioned, the advantages of MCNPX are a complete set of nuclear libraries, very good neutron 
simulations and a card to calculate fission energy deposition. Therefore MCNPX is still the best simulation tool for 
the purpose of the present study. In the absence of detailed designs of nuclear weapons and a complete simulation 
package, the present feasibility study is meant to be an order of magnitude estimation of fission energy deposition 
from annihilation-induced products. 

The geometry of the primary stage of the W87 thermonuclear warhead used in the MCNPX simulations is illustrated 
in Fig. [3l Due to the spherical symmetry of the geometry, a uniform disk source is adequate to simulate random 
radiations on the core. On a 2 GHz Core Duo 2 processor, a particle history of an external takes 2-3 hours to 
complete on one core. Primaries of other particle types may take only 2-3 minutes to finish. In order to be fair in 
the comparisons of all particle types, the same number of particle history (NPS) is used for all of the simulations, 
namely NPS = 100. Table |TT] summarizes the annihilation-induced fission energy deposition on the plutonium core of 
a mocked- up W87 warhead. The unit of fission energy deposition output by MCNPX is MeV/g/p (where p stands for 
"primary") and is converted to MeV/p by multiplying the output with the mass of the plutonium core 1.0015 x lO'' g. 
The energy deposition on the human body is simulated by shooting the primaries from a point source inside a 
rectangular parallelepiped of 200 x 100 x 100 cm^ made with a material of similar chemical composition of the human 
body as shown in Table [TTTl The simulation results of annihilation-induced energy deposition on the human body is 
summarized in Table HVl The unit of energy deposition in MCNPX is MeV/g/p and is converted to Gray per primary 
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(Gy/p) where Gray is defined as J/Kg. The large errors in Tables HTl and [TVl are partly due to varying slant depths 
when the primaries intersect different parts of the targets. 

A typical accelerator such as MINOS can produce neutrinos (93% ly^, 6% and 1% i/g + I'e) from a beamline of 
2 X 10^'^ protons/s. High energy protons interact in a target on the NUMI beamline to produce pions and kaons, 
which decay into muons and muon neutrinos (tt^ —^ ^'^v^). Muons further decay into electrons and neutrons 
— > e^VeV^i). The energy spectrum peaks at approximately 3 GeV with a long high-energy tail extending to 
120 GeV |26|. Assuming a hypothetical situation in that the peak of the neutrino spectrum is tuned to M2/2 for 
the pole, much of energy is still lost to smearing and the production of by-products. An accelerator of the size 
of Fermilab is not practical for neutrino weapon applications because of the problem of maneuverability and the 
prohibitive cost. According to Table [III the typical fission energy deposition is in the 10^ MeV/p range. Assuming a 
very optimistic estimate of 5% efficiency for energy transfer, neutrinos produced by an optimized version of MINOS 
and its antineutrino counterpart will create enough annihilation-induced radiations to convert the plutonium core in 
Figure [3] into a 1-10 kW reactor. Suppose that the chemical explosive around the plutonium core is C4 which is 
91% RDX (cyclotrimethylene-trinitramine) and has a flash point at 234°. Reference [1| estimates the heat capacity 
of ^^^Pu as 6.557 x 10^^ MeV/g-K while Reference [27'] shows a more detailed temperature dependent heat capacity 
profile for both pure and alloyed plutonium as 30-40 J/mol-K at 20-250 °C. As an approximation, this work takes the 
heat capacities of the plutonium alloy and beryllium to be 8 x 10^^ MeV/g-K and 1.14 x 10^^ MeV/g-K respectively. 
A MINOS-induced 1-10 kW heater can detonate a plutonium core in Figure [3] in 100-1000 seconds. Long heating 
time will likely lead to the heating of the area surrounding the plutonium core. Before the temperature in the primary 
stage reaches the flash point of the explosive, the polyethylene filler may melt and the supporting structure may 
expand. If the position of the fissile shifts relative to its surrounding, the effectiveness of the explosive lenses may 
be affected in such a way that a fizzile (not a full scale nuclear explosion) occurs. In case the nuclear weapon is 
moving, it is tactically advantageous to minimize the heating time because of the technical difficulty of tracking. For 
these reasons, the heating time should be limited to approximately 1 s to detonation. The main disadvantage of a 
conventional accelerator is not just the long heating time but also its size and cost. Fortunately recent breakthroughs 
on tabletop accelerators may one-day be the solution for all these problems. Plasma accelerators utilizing the laser- 
driven Wakefield can achieve an electric field gradient up to 270 GeV/m depending on the electron density in the 
plasma. The electron beam produced by Laser Wakefield Acceleration (LWFA) is stable and collimated [2^ . LWFA 
has some potential limitations. A study shows that electric field breakdown occurs at 13.8 ± 0.7 GeV/m in some 
dielectric [29]. In one report, a mono-energetic electron beam up to 1 GeV driven by a 40 TW laser in a 3.3 cm-long 
hydrogen-filled capillary discharge wave guide has been successfully demonstrated [30|- The peak current achieved is 
up to 300 A (2 X 10^^ electrons/s). At the time of the writing of this paper, the highest peak current published by 
another report is 100 pC/10 fs or 10 kA [3l|,[3l|. Results of a positron laser-plasma accelerator is recently pubhshed [3j| 
so that all of the essential components for the construction of vv beams are already in place. Assuming a typical 
femtosecond laser pulse rate of 30 kHz at the current level of technology, LWFA vv beams of 1-10% efficiency for 
energy transfer can detonate the plutonium core in Figure [3] in 50-500 s. 

The radiations produced by the neutrino counter nuclear weapon can also irradiate human beings. The lethal dose 
for a typical human adult is 10-20 Grays. According to Table ITVl the dosage from annihilation-induced radiations is 
in the lO^^** — 10~^^ Gy/p range. Assuming an average dosage of 10"^'^ Gy/p per beam and the same parameters of 
the LWFA vv beams above, a neutrino tactical weapon takes 10-100 s to achieve the lethal dose per person. 

IV. DISCUSSIONS AND CONCLUSION 

An order of magnitude calculation based on a mocked-up design of the W87 warhead and a set of preliminary 
MCNPX simulations shows that radiations from vv ff reactions at the pole is capable of detonating a nuclear 
weapon remotely. The theory of the neutrino weapon capitalizes on the large cross section of the Z'^ pole at oblique 
angles. If the assumptions of relativistic arguments leading up to Eq. (|24p are incorrect and that the width is 
not substantially reduced as claimed, the efficiency of a neutrino weapon will seriously suffer to the point of being 
impractical. Assuming that the theory is correct, the efficiency of the neutrino weapon will largely depend on our 
ability to control the parameters in Eq. pB)) to satisfy the condition = exactly. Since ^ l/[9^ — M'^]'^ 

in Eq. (fT5|) . an agreement between E^, and i?p with AIz to 4 significant figures implies a 10^^ increase in the cross 
section and so on. In reality, there will always be a smear of neutrino energy around the Z'^ pole. Too much smearing 
will decrease the efficiency of the neutrino weapon. The reaction vv — > vi> also decreases the efficiency of the neutrino 
weapon but this effect is predictable and inconsequential. The heat transfer from the plutonium core to the chemical 
explosives can in principle be reduced by thermal insulation or by changing the design to a gun-type trigger mechanism. 
These fixes will most likely reduce the efficiency of a nuclear weapon and may not protect the weapon from being 
detonated remotely. The neutrino counter nuclear weapon can also be used as a tactical weapon. The advantage of a 
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neutrino tactical assault weapon is that it can hit hard-to-reach places such as deep underground concrete bunkers and 
caves in mountainous areas. Since neutrinos can travel across the globe in the speed of light, a neutrino weapon does 
not allow the time for an early warning system. Since neutrinos cannot be shielded, a neutrino weapon is in principle 
non-defensible. A preliminary analysis of the capability of a neutrino weapon given above is based on the experimental 
designs of plasma accelerators. The order of magnitude estimates show that the current level of technology is already 
meeting the minimum hardware requirements. It is foreseeable that high quality and high intensity e^e~ beams will 
be made available by LWFA and similar technologies soon. The remaining technological challenge is to convert the 
e+e" beams into coUimated mono-energetic high intensity vv beams. Accelerator R&D will be the most critical part 
of the development of a neutrino weapon. Experimental verification of the near-singularity of the pole at oblique 
angles is also very important. Since neutrino propagation through the Earth is affected by neutrino oscillation and 
matter effect, basic neutrino physics research is needed to measure the physical paramters accurately. In order to 
obtain maximal collision efficiency, the vD beams must be focused on a small spatial region. Due to the relatively low 
rate of lethal radiation on a human body, a neutrino tactical weapon is a unlikely candidate for a weapon of mass 
destruction. At the end, it is hoped that the neutrino counter nuclear weapon will make nuclear weapons obsolete so 
that the goal of non-proliferation is achieved peacefully. 
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II V 



(a) (6) 

FIG. 1: Tree-level Feynman diagrams of neutrino-antineutrino (yv) annihilation with (a) a charged current and (b) a neutral 
current. The products of vv annihilation is a lepton-antilepton pair (IV) in the charge current case and a fermion-antifermion 
pair (//) in the neutral current case. The fermion / can be a lepton, neutrino or quark. In the case of quark production, 
mesons will be created by hadronization. 
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(a) (b) 

FIG. 2: Graphical representation of incoming 3-momenta pi and pa and outgoing 3-niomenta pa and p4. The point O lies 
on the ellipsoid E in part (a) of the figure. The point X is fixed by the incoming 3-momenta. Part (a) shows a special case in 
which all the 3-rnomonta lie on the same plain. T and T' are the foci of the ellipsoid. The point M lies on the minor axis of 
E.\i M.= T, Q\ = Q-2. Momentum is conserved on E, namely pi -|-p2 = ps +P4. Part (b) illustrates the Cartesian coordinates 
used to define the set of 3-momenta. 
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7.5 cm 



FIG. 3: A mock-up of the primary stage of the W87 thermonuclear warhead with the explosive lenses removed. The core is made 
out of Plutonium alloy (99.4% ^^^Pu and 0.6% natural gallium). The density of plutonium alloy is 15.86 g/cc. The neutron 
reflector is made of beryllium. The density of Be is 1.85 g/cc. The tritium-deuterium booster has a density of 0.00175 g/cc. 
The source is / € {e=^, 7r=^, K"^, Kl, Ks}. 
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TABLE I: Neutral vector and axial vector couplings in the GWS Model '5]. 



Cv Ca 



e~, r~ -i+2sin6'„ -i 

M, c, t I ~ I ^^11 ^™ \ 

d, s, b — i + |sin^u, — i 
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TABLE II: Fission energy deposition Ec from annihilation-induced products used as primaries p on the plutonium core of 
the primary stage of a mocked-up W87 thermonuclear warhead. The geometry is illustrated in Fig. |3] The number of particle 
history is 100. The energy of the primaries is Ef = 45 GeV. 



p 






(10'* MeV/p) 


e± 


2.30 ±0.36 


7r± 


5.80 ±0.59 


if± 


2.59 ±0.39 


Kl 


2.11 ±0.34 


Ks 


1.58 ±0.25 
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TABLE III: Approximate chemical composition of the human body. The average mass and density of a human body is taken 
to be 75 Kg and 1.4 g/cc respectively. 



Element Relative Atomic Composition 

O 65 

C 18 

H 10 

N 3 

Ca 1.5 

P 1.2 

K 0.2 

S 0.2 

CI 0.2 

Na 0.1 

Mg 0.05 
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TABLE IV: Energy deposition Eh from i/v annihilation-induced products used as primaries p on the human body. The chemical 
composition of the human body is given in Table IIIII The number of particle history is 100. The energy of the primaries is 
Ef = 45 GeV. 



p 






(Gy/P) 


e± 


(2.22 ±0.10) X 10"^^ 


7r± 


(5.03 ±0.63) X 10"^'' 


if± 


(4.06 ± 2.07) X 10"^'' 


Kl 


(3.62 ± 1.03) X 10"" 


Ks 


(3.39 ±0.96) X 10-" 



